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Understanding mechanisms for energy dissipation from nanoparticles in 
contact with large samples is a central problem in describing friction mi- 
croscopically. Calculation of the reduced density matrix appears to be the 
most suitable method to study such systems that are coupled to a large en- 
vironment. In this paper, the time evolution of the reduced density matrix 
has been evaluated for an arbitrary system coupled to a heat reservoir. The 
formalism is then applied to study the vibrational relaxation following the 
stick-slip motion of a small adsorbate on a surface. The frequency depen- 
dence of the relaxation time is also determined. 
PACS: 44.10.+i;62.20.Qp;63.22.+m;67.55.Hc 







I. INTRODUCTION 



Friction and its microscopic origins have been intensively investigated in recent years.S 
Progress in atomic force microscopy (AFM) has made possible precise force measurements 
at the atomic scale.! These data, together with simulations, theory and ab initio force cal- 
culations have helped to provide a detailed picture of friction. Much of this work is 
concerned with the nature and rate of energy transfer from lubricant layers or asperities 
into bulk substrates. The understanding of such processes has implications as far-ranging 
as the design of solar collectors,^ where collected energy must be transferred to the bulk 
before reemission to the surroundings; to the understanding of reactions in living cells cat- 
alyzed on small surfaces,!!! as well as biomolecular motorsjl! and to the design of lubricants 
for industrial purposes.0 In many cases, friction is the result of energy dissipation by mul- 
tiphonon processes. We address in this article energy dissipation in a stick-slip model for 
friction involving small and sparsely distributed adsorbates on a surface. We determine the 
nature and calculate the rate of damping of the adsorbed molecule, nanoparticle or asperity 
following the slip step. 

Relaxation processes of adsorbates have been studied both theoretically and experi- 
mentally (see eg.0^0). In most of these works, the authors used either phenomenological 
approximations!!! or assumed that the system under study is at least at quasi-equilibrium 
and used equilibrium theory to study their propertiesBElSS. In this article, a Redfield 
Theory-like approach is developed (for the derivation of the Redfield theory and some of 
its applications see, e.g., Refs. [&"§!]) to calculate the time evolution of the reduced den- 
sity matrix of an adsorbate on a substrate within the framework of a stick-slip model for 
friction. 

Consider an object adsorbed on the surface of a sample, a metal or insulator, with 
vibrational frequency Q. In general there are two possible decay modes: 1) it can create 
electronic excitations in a metal, eg., create electron-hole pairs, or 2) it can create phononic 
excitations. In this article our interest will be in phononic dissipation. If Q ~ uujq where 
ujq is the maximum phonon frequency of the sample, the excitations can decay only by the 
creation of n phonons in the sampled!. For large n, this contribution is in general negligi- 
ble. Thus for small particles on a solid substrate, for which the vibrational frequency of 
the adsorbate might be much higher than the band of substrate mode frequencies, energy 
dissipation from the adsorbate into the substrate may be quite slow. For such cases, an- 
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harmonic coupling between the adsorbate and substrate gives rise to energy flow into the 
substrate. For example, for CO adsorbed on Cu, the Cu-CO stretch vibration is Q ~ 1.5 Uq 
and decay by creating two phonons might be an important mechanism for the vibrational 
damping of the molecule. We will address the problem of phononic energy dissipation in 
friction following stick-slip motion of an adsorbate whose vibrational frequency lies above 
the band of substrate frequencies, like CO on a Cu surface. In Ref. [@] two and three 
phonon contributions to the dissipation of the Cu-CO stretch vibrations was studied using 
the Golden Rule formula. We will study this and similar systems here using a more general 
Redfield theory-like approach for various f2's to understand the dependence of the dissipa- 
tion rate on the coupling between the sample and the adsorbed atom. We assume that the 
adsorbates are sparsely distributed so that we can neglect interactions among them. 

The organization of the paper is as follows: In Sec. II, we calculate the time evolution 
of the reduced density matrix which allows one to take into account all non-equilibrium 
properties of a system and also takes into account possible coherence and incoherence 
effects (for the properties of density matrices see, e.g. Ref. [@]). Possible limitations on 
the obtained evolution are also discussed. In Sec. Ill a model system is proposed which is 
analyzed and solved in Sec. IV. Concluding remarks are given in Sec. V. 

II. EVOLUTION OF THE REDUCED DENSITY MATRIX 

In studying the dynamics of systems coupled to an environment, it is most natural to 
use the Reduced Density Matrix (RDM) formalism. The reduced density matrix, pn, is 
obtained from the density matrix, p, of the system plus the environment by taking its trace 
over the environmental degrees of freedom, 

p R (t) = Tr b p(t) , (1) 

where Tr^ denotes trace over the degrees of freedom of the environment. In the following, the 
index R is omitted whenever there is no possibility of confusion. Once the time evolution 
of the RDM of the system is known, the time evolution of the expectation value of any 
observable, O, of the system can be obtained as: 

(0)(t)=Tr[p R (t)0], (2) 

The time dependence of the RDM of the system can be obtained from the definition of 
the RDM in Eq. (0). Let 
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H = H s + H r + Hi n t , (3) 

where Ti s , H r are the system and reservoir Hamiltonians, respectively, and Hi nt describes 
the interaction between them. We will assume, without loss of generality, that 

n mt = Y,QsF s (4) 



s 



where Q S (F S ) acts only on the system (reservoir) degrees of freedom. The time evolution 
of the components of the RDM is given by 

a'P' 

where hu a p = e a — tp and the tensor R aa ';/3i3'(t) is defined as 

Rocs-wit) = J2P(Ej)(at\S(t)\a'j)((3'3\&(t)\Pi) - <W<W, (6) 

ij 

where the scattering matrix, S(t), is defined as 
S(t) = ei Hot e^ m 

= 1 - ^ J* dt"H int (t') + f dt ij* dt 2 H int {t 1 )H int {t 2 ) + ... (7) 

Here, H is defined to be H = H r +H s . Also H in t(t) = e% nt>t H int e-% nt,t , and = \i)®\j) 
with, 

n r \j) = E J \ 3 ). (8) 

In the following Greek (Latin) letters will denote the system (reservoir) degrees of freedom. 
In deriving this result it is assumed that the reservoir is always in equilibrium so that the 
density matrix of the whole system can be factorized as 

( 1] \p(t)\5k) = 5 3kPRl sP{E j ) (9) 

where the diagonal density matrix elements of the reservoir are defined as 

P{E 3 ) = — , (10) 

where Z = E, e _/3 ^. 



Until this point, the only assumption made is that the density matrix of the whole 
system is factorizable which resulted in a linear equation for the components of the RDM. 
The applicability of this approximation should be studied carefully. This assumption is valid 
only if there exists a weak coupling between the system and the reservoir so that the tensor 
product states \aj) can be considered to be nearly the eigenstates of the whole system. 
If there is strong coupling between the system and the reservoir, or if the "reservoir" is a 
finite one, the density matrix of the whole system in general cannot be factorized and one 
has to do without this simplifying approximation. This approximation should be valid for 
our study of phononic energy decay from an adsorbate coupled anharmonically to a large 
substrate. 

Now, the main task is to find a suitable approximation for the tensor R a a'\pi3'{t). Once it 
is known, the time evolution of the RDM can be calculated. Unfortunately, the expression 
obtained by straightforward application of the second order expansion of the S matrix yields 
a result which is valid only if the time, t, is sufficiently short. To overcome this difficulty we 
used an iterative scheme in which we calculated the initial RDM and then evolved it one 
step in time; taking the evolved RDM as the initial RDM, we evolved it one step further. 
At each step, we evolved for a short enough time. Since energy is not conserved for finite 
times, one has to impose energy conservation by hand. For this reason, the matrix elements 
of 'Hint coupling states of different energies are neglected. The calculations are similar to 
the ones done in scattering theory with the result 

Pa p(t + At) = p aP (t)e-^ At + E Rc*#p>(te)p*p(t)e- iu <#* t , (11) 

a'P' 

Raa';f3(3'(At) = -TW^a-a' fi-fi 1 ( QT Qs'^jss'i^aa') — 
n \ss' 

A ss' 7 

- Uc* E E QTQTjss'm) , (12) 

A ss' 7 / 

j*(u,) = M£p(E J )F*F i J', (13) 

kj 

where the prime on the summation in Eq. ([T3|) indicates that the sum should be carried 
out over states for which hu = E^. 
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III. THE MODEL HAMILTONIAN 



Consider an object adsorbed on a surface. Let M be the mass of a reservoir atom and 
m be the mass of the adsorbed object. Assume that the small adsorbate is bonded to a 
single atom of the sample and the interaction between the sample atom and the adsorbed 
atom is described by the Morse potential, 

U(u -v) = E { e - 2a(u - v) - 2e' a(u ~ v) } , (14) 

where u and v are the vertical displacements of the adsorbed atom and the sample atom, 
respectively. E is the binding energy of the adsorbed atom and a can be related to the 
vibrational frequency, Q, of the adsorbed atom through 

Q =fe) ' (15) 

Expanding the potential and retaining the lowest order terms, we get 

n in t = Auv + Buv 2 (16) 

where 

A = -2E a 2 (17) 
B = -3E a 3 . (18) 

For Q > uq, the uv-term makes no contribution to dissipation since it does not conserve 
energy, and we only have the wy 2 -term. For the other case, Q < u , the contribution to the 
decay rate from the uv term is generally much larger than from the uv 2 term. The decay of 
the vibrational excitation of the adsorbate to the substrate for the case where Q < u using 
only harmonic coupling has been studied by exact diagonalization of the Hamiltonian0. 
The calculated value for the decay rate in Ref. p3| is two orders greater than the value we 
have calculated in Sec. IV, where we used only the uv 2 term. In this article, we consider only 
effects of the uv 2 term. We thus focus on energy dissipation for high-frequency adsorbates, 
where Q > ujq. Then the full phononic Hamiltonian of the system we study becomes 



H = htttfb + + Buv 2 , (19) 



ko- 

where uj^ are the frequencies of the sample phonons with wave vector k and polarization 
vector e a , and b and &k<r are the annihilation operators for the phonons at the adsorbed 
atom and the phonons in the sample, respectively. 
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IV. NUMERICAL ANALYSIS AND DISCUSSION 



In order to construct the initial density matrix, consider the following situation: Assume 
that two samples, Sample I and Sample II, are moving on top of one another with an 
adsorbed layer on the bottom one, and there is no direct interaction between the samples, 
as illustrated in Fig. [I]. Consider the case when the coverage of the adsorbed layer is so 
low that the interactions between the adsorbates can be neglected, in which case one can 
treat each adsorbate independently. During the motion of Sample II, the atom adsorbed 
on Sample I will be pushed and released. If there is a step dislocation on the bottom 
surface of Sample II, the atom will be adiabatically pushed down, due to the wedge shape 
of the surface, displacing it from its equilibrium position and storing energy in it. Then 
it is suddenly released. After its release there is no interaction of the adsorbed atom with 
Sample II. This stick-slip model is relevant for energy dissipation through phonons in dry 
sliding friction or lubrication, as well as for the vibrations of the adsorbed species. The 
character of the contribution of such a mechanism to friction between the bodies would 
depend on the rate of relaxation of this non-equilibrium situation. 

Initially, the density matrix of the system plus reservoir is the equilibrium density matrix 

e -/3(e a +E j ) 

p = E — v — Mto'l. (2°) 

aj 

where Z = J2aj e~^ ( - ta+Ej \ a denotes the number of phonons of the adsorbate and j is a 
multiple index describing the number of phonons in each mode, kcr, of the sample. 

Adiabatically displacing the atom would not cause the atom to go out of equilibrium. 
The density matrix will still be diagonal with the same diagonal elements but in the dis- 
placed basis 

p = E — v — Wi){ctj\ (2i) 

aj 

with the same Z. The displaced harmonic oscillator states, \a'), are defined as 

\a') = e^\a) = J2ca>f3\(3), (22) 

P 

where s is the displacement of the oscillator and p is the momentum operator of the adsor- 
bate. When the adsorbate is suddenly released, the density matrix does not change, but 
now, in the absence of the external force due to Sample II, the density matrix is no longer 
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diagonal in the energy eigenstates, and the adsorbate is out of equilibrium. Denoting the 
RDM of the system right after it is released by p(0 + ), we have 



p(0" 



e ~l3(e a +E 3 ) 

Tr b E — z — I a 'fi I 

aj 

E^(0)l/?)(7l 

Pi 



(23) 



where p j g 7 (0) = Z 1 Y J aC a 'pc* a , e 



-/3e a 



Following Ref. a, we take 



(24) 
(25) 



where Af is the number of the Sample I atoms and is the polarization vector of mode 
kcr. As is pointed out in Ref. [H], this expression for v does not account for the surface 
which might reflect bulk phonons, and also does not take into account any surface phonons. 
With these definitions and choosing 



we obtain 



+ 



(CO 



2M 



F 1 = v\ 
Qi = Bu , 



duj'g(uj')g(uj' - uj)A 



(26) 
(27) 



(n B (oj') + l)n B (u)' - to) 



uj'(u' — uj) 



2M t 
+ 9(-u) 

( n x 

+ 



duj'g(uj')g(uj — uj') 



n B (uj')n B (uj - to') 



uj'{uj — UJ 1 ) 



+ 



2M 



duj'g(uj')g(uj' - uj)A 



(n B (uj') + l)n B (uj' - uj) 



uj'(uj' — uj) 



2M 



, , , , w n B (u )n B {uj - uj 

duo g{uj )g(uj — to )- 



(28) 



uj'(uj — UJ 1 ) 

where the integration region in each integration is where the density of states is nonzero 
and uj is positive. In this result we have assumed the thermodynamic limit and neglected 
O(jj) terms. In this study, g{uo) is represented by the Debye density of states 

, , SuJ 2 



UJ D 
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UJ 



UJ D 



(29) 



where up is the Debye frequency, which is coq, the maximum frequency of the substrate. 
In order to obtain Eq. (2S) from Eq. flT3|), the summations over states are converted to 



integrations over energies and the integration region is chosen so that only a small energy 
violation, Au, is allowed, which is assumed to satisfy AuAt = 1 from the energy-time 
uncertainty relation. If one compares Eq. ( p8| ) with similar results found in the literature 
(eg. Ref. [H]), there is an extra factor of n which arises because of the assumption that At 
is large enough so that one can take the limit t — > oo in certain integrals. This factor is not 
related to the formalism but is just related to the evaluation of Eq. (|T3|). 

The final result can be compared with the results in Ref. [@] . In Ref. |@] it is assumed 
from the beginning that only the diagonal element of the density matrix corresponding to 
the first excited state, pn, is non-zero. In this case, the contribution of the other elements 
of the density matrix can be neglected in the evolution of pn, and we obtain 

p n {t + At) = p n (t) + R n . u (At) Pn , (30) 

which yields a decay rate — R ^L U , which is nothing but the result derived in Ref. |@] using 
the Golden Rule (other than an overall factor of n, as discussed earlier). This feature is 
quite general in the sense that as long as just the first few elements of the density matrix are 
important, and for sufficiently low temperatures, the results obtained using this formalism 
and those obtained by the Golden Rule are almost identical. 

For the numerical data, we have used the following values: huo = 37.6 meV; M = 
28 amu; m = 28 amu; hVt Q = 46 meV; E = 1.8 eV; F = KT 10 N and T = 300 K. Here F 
is the maximum vertical force applied to the adsorbed atom and is related to the vertical 
displacement, s, through 

F 



mVt 2 



(31) 



We vary the adsorbate frequency Q from 0.2Q Q up to 1.52f2 , and take the iteration step 
to be At = fl- 1 . 

In Fig. ||, we have plotted the decay profiles, -|^y, i.e., the time-dependent energy in 
the adsorbate mode over its initial energy, for various values of the adsorbate frequency, 
Q. In each case the range of the time axis corresponds to 300 iterations, each iteration 
corresponding to a time of The exponential character of the decay is clear. For the 
numerical calculations, we used only a finite, 16 x 16, part of the infinite density matrix. 
This caused the matrix elements at the edges to evolve incorrectly. However, as long as 
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they are negligible compared to the matrix elements corresponding to the first few excited 
states, this does not affect the general profile of the time dependence of the energy, which is 
mainly determined by the evolution of the first few diagonal elements of the density matrix. 
In most cases after 300 iterations the matrix elements at the edge become nonnegligible. 
In all cases, we found that the excess energy can be fit almost perfectly to the expression 

AE(t) = AE(Q)e-r , (32) 

where AE(t) = E(t) — E{oo). Here, r is the decay time constant (or the relaxation time). 

In Fig. H, the dependence of the decay time, r, on the adsorbate frequency, Q, is 
shown. In the graph, the frequencies are given in units of Qq, the representative value 
of the CO-Cu frequency given above. We note that our calculations only include cubic 
anharmonic coupling of the adsorbate to the substrate. This is reasonable when Q > U£>. 
When the adsorbate frequency is smaller than the Debye frequency, indicated in the figure, 
the plotted r will be too large, since we have neglected linear coupling terms, which are 
important when Q < oj£>. We nevertheless include the regime where Q < ujb to show how 
the cubic anharmonic contribution to r varies in this regime, though the focus of this study 
is where the adsorbate frequency lies outside the band of substrate frequencies. 

We see that both for large Q and small Q, r diverges. In the large Q limit, the reason 
is due to the phase space factors; the two phonons created or absorbed have to be in a 
band of width 2ujd — ^ which goes to zero as Q — > 2lod- For Q > 2ud, the adsorbed atom 
cannot decay through the emission of two phonons and one has to consider three or more 
phonon processes. In the small Q limit, the coupling constant B becomes very small and 
the system behaves almost as if it is isolated, and cannot decay. In Fig. [| one also sees 
that in the region Q ~ ujd there is an interesting variation in r. For Q < ujd there is a 
contribution to the decay process whereby the adsorbate absorbs a low energy phonon from 
the sample which unites with an adsorbate phonon to emit a high energy phonon into the 
sample. This process is absent when Q > ud, which causes r to rise when Q exceeds ujd- 
We note that this rise, followed by a drop in r with larger Q, is not affected by the neglect 
of linear coupling terms since this interesting nonmonotonic variation of r with Q occurs 
where the adsorbate frequency is higher than the band of substrate frequencies. 
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V. CONCLUDING REMARKS 



In this work, we developed a formalism based on the reduced density matrix to study 
the dissipation of excess energy of an adsorbate into a substrate in the framework of a 
stick-slip model for dry friction. Our focus was on adsorbate frequencies that lie outside 
the band of substrate frequencies. In this case, anharmonic coupling between the adsorbate 
and the substrate controls the rate of energy dissipation, and we introduced cubic terms 
explicitly in our calculations. This work presents an extension and generalization of earlier 
work in which a Golden Rule formula was used to calculate the rate of energy dissipation 
from an adsorbate following the slip step.@ In that work, the adsorbate was linearly coupled 
to the substrate, and anharmonic terms were not considered. In this article, we found that 
excess energy decays essentially exponentially from the adsorbate into the substrate, using 
the reduced density matrix formalism with anharmonic coupling and a range of parameters 
that include those for CO adsorbed on a Cu substrate at room temperature. Previously, 
we assumed exponential decay in adopting the Golden Rule. Here, we have developed a 
more general formalism that yields exponential decay for the systems we have investigated 
at room temperature. 

Our work provides an atomic scale theory explaining how and how fast energy dissipates 
in a stick-slip model for friction when adsorbates are both small and sparsely distributed. 
When the adsorbate is small, the adsorbate frequency may lie outside the band of substrate 
frequencies, and the anharmonic terms that we explicitly take into account here will lead 
to energy dissipation. Because the adsorbates are distributed sparsely over the substrate, 
we can neglect the interactions among adsorbates in addressing the energy dissipation from 
each one. The model we have presented fits very well to friction force microscopy. Our 
work complements a number of important recent studies that address energy dissipation by 
phonons through films of lubricantsilE and through the collective motion of asperities,^ 
systems where the interactions of the physisorbed or chemisorbed adsorbates give rise to 
collective excitations involved in friction. The theory we have presented provides a clear 
picture of phononic energy dissipation through a nanoparticle, and allows quantitative 
analysis on this scale. 
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FIG. 1. An adsorbate between the surfaces of two samples, one of which moves with a velocity 
v. (a) There is no interaction between Sample II and the rest of the system, (b) The adsorbate is 
squeezed, absorbing some of the translational energy of Sample II. (c) The adsorbate is suddenly 
released causing it to oscillate and the interaction between Sample II and the rest of the system 
can again be neglected. 
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